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Evidence for a Dual-Source Mechanism of THz Radiation from Rectangular Mesas of 
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The THz radiation emitted from rectangular mesas of single- crystal Bi2Sr2CaCu2 0s+5 was 
studied using angular distribution measurements and Fourier transform infrared spectroscopy. 
Unlike the recent theoretical predictions, the results provide strong evidence for a dual-source 
mechanism in which the uniform and non-uniform parts of the ac- Josephson current act as 
electric and magnetic current sources, respectively. The latter synchronizes with cavity modes 
of the mesa with integral harmonics of the fundamental radiation. 
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The recent discovery of coherent electromagnetic 
waves at terahertz frequencies (1 THz=10 12 Hz) from the 
intrinsic Josephson junctions (IJJs) within a mesa fabri- 
cated from high-quality single crystals of the high tran- 
sition temperature T c superconductor Bi2Sr2CaCu20s+^ 
(BSCCO), presently denoted as STAR (Stimulated Tera- 
hertz Amplified Radiation) emitter, has triggered a great 
deal of research interest. Besides fundamental studies of 
the physics and chemistry, there has been a great inter- 
est in the development of possible compact and all solid- 
state THz radiation sources, which may be very useful 
for applications in science and engineering fields such as 
medicine, diagnostics, pharmaceutical development, bio- 
science, ultrahigh-speed communication, environmental 
studies, security issues, and various types of nondestruc- 
tive and noninvasive sensing and imaging, etc. 1 "* Some 
basic features of the STAR emitter were previously re- 
ported: 2,3 ^ the emission is stable, continuous and power- 
ful (~ /iW), and its spectral width is less than the res- 
olution limit of the Fourier transform infrared (FTIR) 
spectrometer, 7.5 GHz. By varying the applied dc volt- 
age V across the stack of N IJJs, the frequency of the 
THz radiation appears to be locked by a cavity resonance 
mechanism, which fixes it to that of one of the standing 
wave modes inside the mesa. However, the precise nature 
of the radiation mechanism as well as the fundamental 
features were not understood yet. Moreover, to develop 
the STAR emitter with greater power for useful appli- 
cations, it is first important to understand the physical 
nature of the emission mechanism experimentally. 

Theoretically, a unique radio engineering approach to 
the mechanism was made by Pedersen and Madsen, 4 ) 
whereas models based on numerical simulations of the 
coupled Sine- Gordon equations for stacked Josephson 
junctions were widely used by others. 5 14 ) Among them 
some models 5 10 ) assumed a magnetic current radiation 
source (or cavity mode) only without taking into account 



the electromagnetic fields at near-field, and also neglect- 
ing the substrate effect. Here we provide clear experi- 
mental evidence contradicting those model calculations. 

We measured the angular dependence of the far-field 
radiation from various rectangular mesas. Our experi- 
mental results are inconsistent with the predictions of 
either an electric or a magnetic surface current sources 
acting alone. Instead, they strongly suggest a dual-source 
THz radiation mechanism: 15 ) the uniform part of the 
ac- Josephson current acts as an electric surface current 
source, and the inhomogeneous part of the ac- Josephson 
current sets up a displacement current that excites a cav- 
ity resonance mode in the mesa, which locks the radia- 
tion frequency and acts as a magnetic surface current 
source. By adjusting the relative amplitude and phase of 
the two source currents and accounting for the substrate, 
excellent agreement with experiment is obtained. Fur- 
thermore, from FTIR analysis of the dc current- volt age 
(I-V) characteristic, the locking onto the resonant cavity 
mode appears to occur by the sequential synchronization 
of each of the IJJs. 

High-quality single crystals of BSCCO were grown by 
the conventional traveling solvent floating zone method 
using a modified infrared image furnace. 16, 17 ) A piece of 
crystal was annealed overnight to an appropriate dop- 
ing level at 600°C in Ar+0.1% O2 atmosphere. The re- 
sulting crystals are underdoped with T c ~ 75 - 86 K. 
Two types of rectangular mesas were prepared using ei- 
ther conventional photolithography or focused ion beam 
milling technique. Photographs of a cleaved piece of a 
pristine single crystal, an atomic force microscope (AFM ) 
image of a rectangular mesa after processing and the fi- 
nal form of the mesa are shown in Figs. 1(a), 1(b) and 
1(c), respectively. From the AFM measurements shown 
in Figs. 1(b) and 1(d), its length was L ~ 400 /im, width 
w = 77 A ± 4.5 jam and thickness d — 1.2 /im. 

The temperature dependence of the normal state c- 
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Fig. 1. (a): a photographs of a piece of single crystal before pro- 
cessing, (b): an AFM image of the mesa used in this experiment 
made by photolithography and argon ion milling technique, (c): 
the final form of the mesa used here, and (d): a view of the mesa 
cross-section along the A-A line. The top and bottom widths are 
w = 72.9 /im, and 81.9 /im, respectively. The top of the mesa is 
covered by a thin Au electrode layer. The Au electric lead wire 
with 10 /im in diameter is attached by silver paste at the upper 
end of the gold layer shown in (c). 

axis resistance Rj(T) and the I-V characteristics of 
the mesas were measured as described previously. 2,3 ^ 
The frequency spectra of the emission were measured 
by a FTIR spectrometer (JASCO FARIS-1) with a Si- 
bolometer detector. The angular dependence of the ra- 
diation was measured by rotating the mesa sample in- 
crementally relative to the detector. The solid angle of 
the radiation from the mesa to the detector is adjusted 
to have an angular resolution better than ± 2.5°. The 
sample is cooled by a liquid 4 He flow cryostat and the 
temperature is controlled to within ± 0.2 K. 

The angle 9 dependence of the radiation intensity X 
in the xz-plane, or <j) = 90°, from a similar rectangular 
mesa with L ~ 400 /im and width w ~ 60 /im are 
presented in Fig. 2 (a), where 9 and <ft are defined in Fig. 
2(b). As clearly seen in Fig. 2(a), the observed radiation 
is very anisotropic. When L/w ^> 1, as for this mesa, 
T{9, <p) is strongest in the xz-plane. The data shown in 
Fig. 2(a) are typical of data obtained in the xz-plane 
from several mesas. Similar results in the xy-plane of the 
mesa are shown in Fig. 3(a), where the X(#, 0°) data are 
mostly several times weaker and exhibit greater scatter 
than the 1(0,90°) data. 

The general observations are summarized as follows: 
the maximum intensity I max occurs in the xz-plane 
(</> = 90°) at 9 - ±30° from the mesa top (0 = 
0°), where X(0, <fi) is a local minimum of 1(9, 0), with 
X(0, 90°)/2 max « 0.4 - 0.7, depending on the mesa mea- 
sured. This rather large sample-to-sample variation in 
X(0, 90°)/2 max is not due to experimental error, but 
arises from specific differences in the properties of the in- 
dividual mesas. For example, the primary radiation may 
occur at slightly different regions in the mesa. The thin 
Au electrode layer may affect the radiation. In addition, 
the spatial inhomogeneities in the ac- Josephson current 




Fig. 2. (a): a polar plot of Z(<9, 90°) normalized at Z(0°, 90°) mea- 
sured in the xz-plane from a rectangular mesa at T=32.7 K. Each 
data point was measured point by point as follows: first, cycle 
the I-V curve enough slowly for only one direction to find out 
the radiation condition, then, stop cycling at the maximum point 
of radiation. The typical values for / and V are 19.97 db 0.02 mA 
and 0.9214 ± 0.0007 V, respectively. The solid line presents the 
respective <j> = 90° portion of the least-squares fits of the dual- 
source model to the data shown here combining the contributions 
from the uniform ac- Josephson current and the inhomogeneous 
displacement current exciting the (001) cavity mode, with cor- 
rections for a superconducting substrate using Model I. 15 ) The 
dashed curve presents the corresponding cavity model fit for a 
sample suspended in vacuum (see text), (b): a sketch of the co- 
ordinate system. 18 ^ 

due to sample heating and stoichiometry variations, etc., 
which are crucial for the cavity mode excitations, 15 ) are 
likely to be sample-dependent. The specific trapezoidal 
cross-sectional shape of each mesa as seen in Figs. 1(b) 
and 1(d) may also affect its angular dependence. How- 
ever, the precise reasons for this variation are not yet 
well understood. 

If the radiation were simply induced by the fundamen- 
tal cavity resonance mode, as expected from capacitor 
patch antenna theory 18, 19 ) and widely predicted for rect- 
angular BSCCO mesas, 5 10 ) 1(9, <j)) would be maximal at 
9 = 0°. If a uniform ac- Josephson current were the pri- 
mary radiation source, X would vanish at 9 — 0° and 
be maximal near to 9 = 90°, as for simple dipole radia- 
tion. 18 ) Apparently, the experimental results contradict 
both simple explanations. 

Second, the observed 1(9,90°) diminishes strongly as 
9 — > 90° (parallel to the afr-plane). This fact contra- 
dicts the simple uniform ac- Josephson current model, 
for which a maximal intensity is expected at or near to 
9 — 90°. If the refractive index, n, were unity for the 
mesa as for the vacuum exterior, the radiation from a 
cavity mode would vanish at 9 = 90°, as observed. How- 
ever, for the appropriate value n = 4.2, the predicted 
X(90°,90°) from the simple cavity resonance model is 
~ 87% of its maximum value, 15 ) as shown in Fig. 2(a), 
unlike the experiments. 

Third, radiation, with an order of magnitude weaker 
than the fundamental, at several integral higher harmon- 
ics of the fundamental frequency is generally observed, 
as the previous report, 3 ) but subharmonics were never 
observed. Fourth, small but clear minor lobes are also 
observed at 9 ~ ± 75° as seen in Figs. 2(a) and 3(a), 
although their integrated intensity appears to be only a 
few % or less of the total. These lobes may arise from 
higher harmonics, but the spectrum of the entire 1(9, <$) 
has not yet been measured. 
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Fig. 3. (a): a polar plot of the observed Z(0,O°) normalized at 
Z(0°,0°) in the xy-plane of the rectangular mesa at T=30.0 K. 
Each data point was measured in a similar manner as that shown 
in Fig. 2(a). The typical values for / and V in this particular case 
are 14.90 ± 0.02 mA and 0.7995 ± 0.0005 V, respectively. The 
solid curve presents the 4> = 0° portion of the overall dual-source 
fit using Model I. 15 ) The dashed curve is the corresponding cav- 
ity model fit. (b): the 3D plot of the predicted T(0,(f>) from the 
least-squares fit. 

These results can naturally be understood by combin- 
ing both radiation sources, since there actually exist both 
the displacement current and the real ac- Josephson cur- 
rent flowing across the mesa, which have both spatially 
uniform and non-uniform parts. Problems in fitting the 
data near to 6 = ±90° can be overcome by an appropri- 
ate model of the superconducting substrate. 15 ) 

It was shown previously that the fundamental fre- 
quency v satisfies the ac- Josephson relation, 2 ' 3 ) v = 
j^-V/N, where e is the electronic charge and h is the 
Planck's constant. Since w is comparable to the wave- 
length of the plasma waves in the mesa, cavity reso- 
nances may occur as standing waves. Empirically, the 
wavelength A = 2w of the lowest energy state in the 
mesa, leading to the condition v = co/nX = cof2nw, 
where Co is the speed of light in vacuum, n = \fc, and e 
is the dielectric constant of the junctions. This empirical 
relation appears to work very well in many samples with 
different widths and L/w ratios, as shown in Fig. 4. From 
the slope of the line in Fig. 4, n = 4.2 is obtained fairly 
accurately, corresponding to the THz-frequency dielec- 
tric constant e = 17.6, which is about 50% larger than 
the value (e = 12) obtained by infrared spectroscopy. 20 ) 
No anomaly in e for frequencies up to 0.9 THz has been 
observed. 

For d <C L,w, the lowest mode frequencies satisfy 
vomp = \/(m/L) 2 + (p/w) 2 , where m and p are inte- 
gers. 18 ) Hence, the lowest frequency mode for L > w is 
expected to be the (010) mode with ^qio = 2 L°y/e ' ^ s seen 
in Fig. 4, the experimental results are definitely consis- 
tent with the (001) mode, not the (010) mode, contra- 
dicting the simple cavity resonance model. This may be 
related to an instability that could arise from the energy 
loss by the penetration of the magnetic field in the z- 
direction into the mesa, since A ~ A c , where A c is the 
superconducting penetration depth in the c-axis direc- 
tion, resulting in higher inductive energy states for the 
(OraO) modes. But whatever the actual reasons for the 
apparent non-resonant (OraO) modes are, the angular de- 
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Fig. 4. A plot of the fundamental frequency versus l/w. The data 
were obtained from samples with 40 /im^ w ^ 100 /xm and 
300 fim^ L ^ 400/im and variable L/w ratios, prepared using 
four different methods from nine different single crystals. Type 
A represents the data from the reference. 2 ) Type B, C and D 
were prepared by FIB, ion milling using metal mask, and pho- 
tolithography techniques, respectively. 

pendence of the far-field radiation is calculated by com- 
bining the output from the two sources, one of which 
excites the (001) mode taking into account the super- 
conducting substrate. 15 ^ Least-squares fits to the data 
shown in Figs. 2(a) and 3(a) were performed using Model 
I, which averages the two source contributions symmet- 
rically about — 0° , while preserving the Love magnetic 
equivalence principle boundary condition n x H = 
on the mesa edge, 15,18 ) and using Model II, which re- 
laxes the boundary condition while preserving A. 15 ) The 
best fit with standard deviation a = 0.122 was found 
for Model I with mixing parameter ai(0) = 0.310, corre- 
sponding to 24% of the overall intensity arising from the 
magnetic (cavity) source. The cross-sections of the best 
dual- and single-source fits are shown by the solid- and 
dashed- curves in Figs. 2(a) and 3(a), and the full radi- 
ation pattern resulting from the dual-source best fit is 
shown in Fig. 3(b). Preliminary but similar results were 
obtained from numerical simulations. 11, 12, 14 ) 

Figure 5 presents the detailed dc I-V characteristics 
and the corresponding output of the Si-bolometer de- 
tector as functions of V, when the scan is made very 
slowly (2.8 hours for a one-cycle measurement). In this 
case many radiation peaks are observed. Examining them 
in detail, one finds an interesting feature. When / is re- 
duced slowly on the return I-V branch, the radiation 
energy starts to build up gradually until the growth is 
interrupted by a jump to a different I-V branch. This 
is more clearly seen in the inset of Fig. 5, in which the 
sawtooth-like radiation peaks are seen over the range of 
V from -0.67 to -0.57 V. From the detailed FTIR mea- 
surements performed at various points in the sawtooth 
region, it turns out that the observed frequency shifts to 
lower frequency as the power increases (i.e., I decreases). 
This can be interpreted as the gradual decrease of the 
ac- Josephson frequency as I decreases, until the mesa 
eventually enters a cavity resonance state. Meanwhile, 
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Fig. 5. The I-V characteristic (connected data-point representa- 
tion) and the radiation power (vertical line representation) de- 
tected by the Si-bolometer at T=30.0 K. The arrows indicate 
the sweeping direction of the I-V curve. The detector was set 
at the angle #=45° and <j)=90° for this particular measurement. 
The inset is a magnified view of the encircled region. 

the radiation energy increases in the mesa, and when the 
power reaches a threshold, / jumps to a different I-V 
branch. Therefore, the mechanism of the strong coherent 
continuous radiation appears to be due to gradual devel- 
opment of the synchronization in the individual or groups 
of Josephson junctions onto the cavity resonance mode, 
rather than a sudden catastrophic transition. This fact 
strongly suggests that the stronger output power may be 
achieved by preventing the I-V curve from jumping. 

The combination of two conditions necessary 
for the resonant emission may be written as 
V^ t =(c /2nw)(l/K J )(2d/c), where V^ 1 is the required 
dc voltage to meet the cavity resonance oscillation, the 
Josephson constant iO=483.5979 (GHz/mV), and c « 
30.65 A is the c-axis lattice constant of BSCCO. This 
reduces to V^ t= 48. 2 (d/w), which surprisingly depends 
only upon d/w. 21 ^ Since an emitting junction must be 
in its resistive state, its resistance Rj and dc current / 
must satisfy Rjl > V^ f . Because Rj(T) has a steep 
negative temperature coefficient, and BSCCO each 
junction is subject to severe Joule heating, the whole 
mesa may thermally be stabilized. On the other hand, 
because of large reduction of Rj the above condition for 
synchronized emission often may not be satisfied, espe- 
cially for over doped samples, due to their considerably 
smaller Rj values. This may explain partially why we 
initially had many unsuccessful mesas. 

In conclusion, we measured the angular dependence of 
the far-field THz emission from mesas of single-crystal 
Bi2Sr2CaCu208+5. Our results provide strong evidence 
for a dual-source mechanism in which the uniform and 
inhomogeneous parts of the ac- Josephson current respec- 
tively act as an electric surface current source and set up 



a displacement current that excites a mesa cavity res- 
onance mode which locks the radiation frequency, and 
acts as a magnetic surface current source. By adjusting 
the relative amplitude and phase of the two source cur- 
rents and accounting for the substrate, excellent agree- 
ment with experimental results is obtained. From FTIR 
analysis of the I-V characteristic, the radiation appears 
to build-up from gradual synchronization of the individ- 
ual or groups of Josephson junctions to the cavity reso- 
nance mode. 
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